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1 Introduction 



Finite temperature simulations of lattice gauge theory suggest a phase 
transition at high energy density and temperature to a new phase of QCD 
matter — the quark-gluon plasma (QGP). In many models, the tempera- 
tures and energy densities reached in central nucleus-nucleus collisions at the 
Brookhaven Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC) at CERN are predicted to be sufficient for QGP formation. 
Demonstrating that such a QGP has formed and understanding its nature 
will be a very important component of the experimental programs at both 
machines. In this paper, we focus on using measurements of the px spectra 
of the Upsilon resonance family to accomplish these tasks, emphasizing the 
possible advantages of the higher LHC energy as compared to RHIC energy. 

The details of the QCD phase transition to the QGP, as obtained in finite 
temperature simulations of lattice gauge theory, are strongly dependent on 
the number of quark flavors, rif, and the number of colors, N. The critical 
temperature, T c , and the order of the phase transition are particularly sensi- 
tive to the values of n/ and N used in the simulations In a pure SU(iV) 
gauge theory, rif = 0, the phase transition is second order, i.e. continuous, 
for SU(2) and first order for SU(3). For the SU(iV) theory, the critical tem- 
perature, T c , is 260 MeV. When dynamical light quarks are included in the 
lattice simulation, the critical temperature is substantially lower and the or- 
der of the phase transition can be altered: for rif = 2, T c pa 150 MeV and the 
phase transition appears to be continuous; for simulations with more light 
flavors, rif > 3, T c is unaltered but the transition again becomes first order. 
This latter conclusion depends on the relative quark masses || . The critical 
energy density at T c is only weakly dependent on the color and flavor degrees 
of freedom; e c ~ 1 — 2GeV/fm 3 is obtained both with and without quark 
degrees of freedom for N = 2 or 3. For a QGP to be formed in an ultrarel- 
ativistic heavy-ion collision, the initial temperature, T , and energy density, 
e(to), of the quark-gluon system created by the collision must be larger than 
T c and e c . 

Part of the LHC experimental program will be devoted to heavy ion colli- 
sions such as Pb+Pb at \/s = 5.5 TeV. At this high energy, the highest avail- 
able for heavy-ion collisions, perturbative QCD processes will be the domi- 
nant factor in determining the initial state of the post-collision quark/gluon 
system. In particular, at early times, t ~ 1/pr < l/Po> semihard produc- 
tion of minijetsQ could be dominant in establishing the initial conditions for 
further evolution of the system [Q . The appropriate value of po is uncertain 
but should be of order p ~ 2 GeV, in which case to ~ 0.1 fm. The crucial 

1 Minijets are jets with pr > po ~ 1 — 3 GeV ||, usually not observable as individual 
jets below p T ~ 5 GeV 0. 
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quantity in determining whether a quark-gluon plasma forms is the initial 
energy density, given by 



e(*o) = ^ , with E T = T AA (0)a jct (T r > . (1) 

Here, R is of order the nuclear radius, Ra, Taa(O) is the nuclear overlap for a 
central collision?] and a^(Ex) is the first pt moment of the minijet production 
cross section in the central unity of rapidity in pp collisions (if we neglect 
nuclear shadowing) using the assumed value of Pq. The initial temperature 
T is proportional to [e(t )] 1//4 - For given p , a- ]et (E T ) is computable from 
QCD once the parton distribution functions (PDF's) for the colliding protons 
and the types of jets participating in the plasma are both specified. If we 
take Pq = 2 GeV, employ MRS D— ' Q parton densities and include only 
gluon jets (i.e. assume a pure gluon plasma) with rapidity \y\ < 0.5, then 
the semihard production of minijets leads to a large initial energy density, 
e(io) — 1170 GeV/fm 3 , and initial temperature, T = 1.14 GeV, in central 
Pb+Pb collisions at \fs = 5.5 TeV, to be compared to e(t ) = 27 GeV/fm 3 
and T = 0.46 GeV at the lower RHIC energy. See 0, i for details. 

The sensitivity of the critical parameters to the choice of PDF set is 
illustrated in Table [l], following 0, §|. MRS D— ' provides the most optimistic 
initial temperature, with the MRS H distributions M (designed to fit early 
HERA data) not far behind. The Duke-Owens-I (DOl) |K|, and MRS DO' 
H distributions that grow more slowly as x decreases yield smaller initial 
energy density. More recent HERA data continues to favor distributions 



(such as the MRS A parameterization [II]) that at small x are similar to the 



H parameterization, i. e. that do not grow quite so rapidly at small-x as the 
D— ' distributions but are certainly much more singular than either DOl or 
DO'. All the numbers quoted in the table will become smaller when reductions 
in parton density due to shadowing in the nucleus are incorporated. 

The sensitivity of the results to po depends upon the PDFs. Increasing 
p from 2 to 4 GeV decreases a^ et {E T ) by factors of 5, 3, 2, and 1.5 for D— ' 
§, H §, DO' [§J, and DOl [TJ distributions, respectively, at ^ = 5.5 TeV. 
However, e(to), eq. (H), also proportional to po, changes by a factor of two or 
less, even increasing in the DOl case. The sensitivity of e(to) to po at RHIC 
is much greater, declining by as much as a factor of 5 between po = 2 GeV 
and po — 4 GeV. Inclusion of quark as well as gluon jets would increase e(t ) 
by roughly 30% at -^i = 5.5 TeV, but if the plasma is a quark-gluon plasma 
with 2 to 3 active flavors, the equation of state changes so that To actually 
declines. 



2 The nuclear overlap function for a collision of projectile A with target B is Tab(^>) — 
J cPsdzdz' pa(s, z)ps{h — s, z') where the impact parameter b is zero in the most central 
collisions. 
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PDF 


#jets 


# gluon jets 


E T (GeV) 


T (GeV) 


dN^ ons only /dy 


D01 


776 


640 


3200 


0.74 


914 


DO' 


1510 


1250 


5440 


0.81 


1350 


H 


3250 


2710 


10300 


0.96 


2130 


D-' 


5980 


5220 


17600 


1.14 


3360 



Table 1: Dependence of initial temperature and particle density on distribu- 
tion function choice at ^/s = 5.5 TeV. The results are taken from Ref. f|. We 
employ TpbPb(O) = 32 mb~ . The distributions employed are Duke-Owens 
I, DOl 0, the MRS DO' and D-' from [|, and the MRS H distributions 
from ||. 



Regardless of the parameters employed, it is apparent that minijet pro- 
duction yields substantially larger initial temperatures than previously pre- 
dicted (for example, in Ref. and thus enhances the possibility of QGP 
production in thermal equilibrium. Only the extent of this enhancement is 
uncertain. For large enough p and/or extensive nuclear shadowing, the en- 
ergy density and temperature might be reduced sufficiently that the plasma 
would not be fully equilibrated, although a QGP might still be formed. In 
any case, it is important to examine the proposed signatures of plasma forma- 
tion for the increased initial energy density and temperature obtained after 
including the effects of minijet production. 

One of the proposed signatures of the QCD phase transition is the sup- 



pression of quarkonium production, particularly of the J/ip |T3|, I14J . In fact, 
J /if) and if)' suppression have been observed in nucleus-nucleus collisions at 
the CERN SPS with ^ = 19.4 GeV @§. In a QGP, the suppression 



occurs due to the shielding of the cc binding potential by color screening, 
leading to the breakup of the resonance. The cc (J /if), if)', • • •) and bb (T, 
T', • • •) resonances have smaller radii than light-quark hadrons and there- 
fore higher temperatures are needed to disassociate these quarkonium states. 
At current energies, the situation for the J /if) is rather ambiguous because 
the bound state can also break up through interactions with nucleons and 
comoving hadrons — QGP production is not a unique explanation of J/if> sup- 
pression [|I7|, [L8|, 0, ^DJ. At LHC energies there will continue to be ambiguity 
in this regard, and, in addition, J/ if) and if)' detection efficiencies may not be 
large PT . 



In contrast, the T resonances are sufficiently tightly bound that it is more 
difficult to break them up as a consequence of ordinary multiple interactions. 
In addition, they can be detected in the decay mode with relatively 

high efficiency at the LHC {e.g. 33% vs. 6% for J /if) in the CMS detector 
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nj). Further, this detection efficiency derives entirely from the central region 
(| 7/ 1 < 2.4 in the CMS detector [pip, where the predictions for the relative 
suppression of the members of the T family in the presence of a gluon plasma 
can be most reliably tested (given that the baryon density is predicted to 
be low and should not play a significant role). However, it was previously 



assumed that the T would not be suppressed by QGP production |TJ[ |22 
since the T is much smaller than the cc and other bb resonances (such as the 
T' and Xb states), implying that a much higher temperature than expected 



from earlier estimates of T is needed to dissociate the T In view of the 
high initial temperature of a gluon plasma, T ~ 1 GeV, when minijets are 
included, a reexamination of this assumption is in order. If T suppression by 
a minijet plasma is predicted, the T may provide a cleaner plasma signature 
than the J / ip because of its weaker interactions with matter. In high energy 
heavy-ion collisions, the T rate should be large enough to test the predictions 
given here and help clarify the initial state of the system. 

The paper is organized as follows. In the next section, we review the 
model of quarkonium suppression in a QGP and introduce two different sce- 
narios for quarkonium breakup at finite temperature. Section 3 presents the 
quarkonium production model we adopt and includes a discussion of the ex- 
pected T production rate. We then provide a detailed examination of T 
family suppression in a minijet plasma as a function of pt and discuss the 
sensitivity of the results to the initial conditions. In section 5, we elucidate 
the steps necessary to extract information on the state of the system. We 
particularly focus on the T'/T ratio and discuss the effect of cascade de- 
cays from higher states on the signal. We also propose using the Z boson 
as a reference process against which T suppression can be tested. In sum- 
mary, we show that the T family may provide an extremely powerful tool for 
diagnosing the existence and nature of the quark-gluon plasma. 



2 Quarkonia in a quark-gluon plasma 



We follow the arguments of Ref. [14"] to discuss the finite temperature 
properties of quarkonium. The interquark potential for nonrelativistically 
bound quarkonium at zero temperature is 



V{r, 0) 



ar — 



(2) 



where r is the separation between the Q and Q. The parameters are a 
GeV 2 



a 



0.471, m c = 1.32 GeV and m b = 4.746 GeV p 



0.192 
This potential 



is modified at finite temperatures due to color screening, becoming 



V(r,T) 



a 



a 



- M (T)r 



(3) 
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The screening mass, fi{T), is an increasing function of temperature. When 
n(T = 0) — > 0, eq. (@) is recovered. At finite temperature, when r — > 0, 
the 1/r behavior is dominant while as r — ► oo one finds that the range of 
the potential decreases with n(T), making the binding less effective at finite 
temperature. Semiclassically the energy is 

E(r, T) = 2m Q + + V(r, T) , (4) 

where (p 2 )(r 2 ) = c ~ 0(1) ||14|| . Minimizing E(r,T) gives the radius of the 
bound state at each T. For /x(T) above a certain critical value, there 
is no longer a minimum and the screening has become strong enough for 
binding to be impossible and the resonance no longer forms in the plasma. 
The temperature Tp at which this happens is that for which fi(Tp) = A*d- 

The values of fip for the bottonium states are given in Table ^|, along 
with the bottonium mass at \i = fip, M(fip). Also shown are the results 
for the charmonium states. We note that \ip and M(fip) are independent of 
T D , but that Tp depends upon the functional form of fi(T). Table ^| shows 
that the T mass is actually larger at breakup than at T = because the 1/r 
contribution to the potential is dominant, whereas the T' and XbO-P) masses 
are smaller due to the fact that the linear term is more important for these 
resonances. 





T 




Xb(lP) 


J/V 




Xc(LP) 


M (GeV) 


9.445 


10.004 


9.897 


3.07 


3.698 


3.5 


r (fm) 


0.226 


0.509 


0.408 


0.453 


0.875 


0.696 


t f (fm) 


0.76 


1.9 


2.6 


0.89 


1.5 


2.0 


M{ji D ) (GeV) 


9.615 


9.778 


9.829 


2.915 


3.177 


3.198 


fiD (GeV) 


1.565 


0.671 


0.558 


0.699 


0.357 


0.342 



Table 2: Properties of the quarkonium states both at zero temperature and 
at T = Tjj (the breakup temperature for each bound state), taken from Ref. 



T4J . The masses at T = are from the solution to the Schrodinger equation 
with the potential of eq. (fj). Note that r is the distance between the Q and 
Q and Tp is the formation time of the bound state. M(/1d) and \ip> are the 
mass at breakup and the screening mass at Tp,] both are independent of the 
actual value of Tp as determined by the functional form of fi(T). 



It is necessary to know the dependence of the screening mass on tem- 
perature before the breakup temperature Tp can be determined. Since the 
behavior of n(T) is unknown for T ^> T c , we consider two scenarios which 
may be expected to bracket the realistic situation. 
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The first is a parameterization based on SU(iV) lattice simulations [23], 

(5) 

which results in the lowest values of To- Although this parameter- 
ization of /x(T) is independent of T c , we will present our results as 
a function of the ratio T/T c taking T c = 260 MeV, consistent with 
SU(3) lattice calculations. For this value of T c , the T breaks up at 



Td/T c fa 1.5. In Ref. [p2], this parameterization was used with T c = 150 



MeV, leading to the estimate To/T c ~ 2.6 for the T. 
The second is an estimate from perturbation theory |24| . 



where the temperature-dependent running coupling constant is 

2 (T\ = 24^ 

9 \TJ (33 -2n / )ln[(19T c /A M s) (T/T c )} ' 1 } 

In SU(3) gauge theory, T c /A^g = 1.78 ± 0.03 |H| . In this case, T D can 
be quite high and depends strongly on n/ and T c . We present results 
using eq. (H) with n/ = 3 and T c = 150 MeV. 

We shall refer to these two distinct choices for fi(T) as the SU(A^) and 3- 
Flavor QGP's, respectively. The temperature dependence for both estimates 
of n(T) is shown in Fig. [I| with he, indicated for T, T' and Xb- The values 
of To are given in Table |3| for both cases. Even in the first case, the early 
estimates of initial temperatures of order T ~ 260 MeV for the LHC im- 
plied that T suppression would be unlikely [Q, and that even T' and Xb 
suppression would not be certain. 

In contrast, the higher 'minijet' value of To ~ 1.14 GeV, combined with 
the results of Fig. [I] and the To values given in Table ^| indicate that, when 
the SU(A^) parameterization is used, the T may easily be suppressed in the 
minijet gluon plasma at LHC energies. However, the 3-Flavor estimate leads 
to much higher breakup temperatures for the Upsilon states and the distinct 
possibility that the T will not be suppressed^ even though the T' and Xb 
could be suppressed. In either case, the suppression pattern depends on the 
space-time evolution of the plasma. Depending upon circumstances, this de- 
pendence could be either an impediment to the extraction of the physics of 

3 Note that T production from T' and Xfc decays will be suppressed if the T' and Xfc 
states are suppressed. 
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2.0 
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SU(N), eq. (4) 

quarks 4- gluons, eq. (5) 



2 4 6 8 

T/T c 

Figure 1: The screening mass as a function of temperature for the scenarios 
of eqs. dH) (solid) and (||) (dashed). The values of /jd at breakup for the 
T, Y', and Xb are also indicated by the dotted lines. 



the QGP or a very powerful tool. In most scenarios, much can be learned 
about the QGP by simply examining the ratio [22| of production rates for 
higher bb resonances to the production rate for the T. However, for a quanti- 
tative determination of QGP suppression on an absolute scale, a benchmark 
such as Z production, to which T production cross sections may be com- 
pared, is needed. 



3 Cross sections, rates, and the color evapo- 
ration model 

Before turning to the details of T suppression, let us determine whether 
or not the T production rate is large enough for any suppression to be observ- 
able. An estimate of the T yield is not straightforward since the formation of 
bound states from the predominantly color octet QQ pairs is nonperturbative 
and therefore model-dependent. 

The color evaporation model p5f has often been applied to the production 
of quarkonium states below the open charm/bottom thresholds. The QQ pair 
neutralizes its color by interaction with the collision-induced color field — 
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T D (MeV) 




SU(iV) only, eq. (§) 


gluons and quarks, eq. @ 


T 


391 


1125 




260 


447 


Xb 


260 


368 




260 


459 




260 


211.5 


Xc 


260 


204 



Table 3: The values of Tp (in MeV) for the two choices of fi(T), eq. (||) and 
eq. (§. We employ T c = 260 MeV for the SU(iV) case and T c = 150 MeV 
in the case of eq. (||). The lower T c used in the 3-Flavor case results in the 
lower Td values for the ip' and \ c . 



"color evaporation". The Q and the Q either combine with light quarks to 
produce heavy-flavored hadrons or bind with each other in a quarkonium 
state. The additional energy needed to produce heavy-flavored hadrons is 
obtained nonperturbatively from the color field in the interaction region. 
Depending on itlq, the yield of all quarkonium states may be only a small 
fraction of the total QQ cross section below the vl = 2mn threshold, where 
mu is the mass of the lightest heavy-flavored meson. The total heavy quark 
cross section below this HH threshold, a, is calculated at leading order by 
integrating over the QQ pair mass from 2m<g to 2m h, 



o is) 



E 



Ami 



linf 



ds I dxidx 2 fi/p{xx) fj/ P (x 2 ) o-y(s) 5(s - x^s) , (8) 



where ij = qq or gg and 0^ (s) is the ij — > QQ subprocess cross section. The 
color evaporation model was taken to next-to-leading order (NLO) using the 
exclusive QQ hadroproduction calculation |26j to obtain the energy, xp, and 
Pr-dependence of quarkonium production J27], ^8 . 



The division of a into heavy flavored hadrons vs. quarkonium and the 
relative quarkonium production rates are both a matter of assumption in 
the color evaporation model. For the model to have any predictive power, 
the relative quarkonium production rates should be independent of projectile, 
target, and energy. Experiment suggests that this is true for the charmonium 
ratios Xc/ ip and ip'/ ip over a broad energy range [^. |3(]. |3l], | 
bottonium data follows this trend: T'/T = 0.53 ± 0.13 and f"/T 
0.06 for proton beams at 400 |£| and 800 GeV [|4|, |35[. These results are 



The available 
0.17 ± 



consistent with the production of T states in pp collisions at the Tevatron 
The color evaporation model also reproduces the energy dependence 
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of charmonium and bottonium production over a wide range of energy as 
well as most of the xp dependence of the charmonium statesQ, provided the 
fraction of quarkonium vs. HH states is assumed to be independent of all 
kinematic variables. 

With these assumptions and our knowledge of the T'/T and T"/T ratios 
at low energy, the normalization of each quarkonium state can be fixed em- 
pirically from data, allowing predictions of the production cross sections at 
LHC energies. The procedure we have followed is described below. 

• As a first step, we assess the accuracy of the model for describing 
existing pp/pp data. Fixed target T data has been generally given 
as the sum of T, T', and T" production. From this data, one can 
determine B(da/dy) y= o where B is an effective dilepton branching ratio 
reflecting the production yield through the T, T', and T" resonances. 



As shown in Fig. ||, taken from Ref. [p7|| , a good fit to the data, up to 



and including ISR energies [B3], B3, |35], |38|, p9f, is obtained with 



where da/dy is computed using the MRS D— ' |J or GRV HO |L0f 



PDFs, taking mj = 4.75 GeV and the renormalization and factorization 
scales set to fi = m T ^ = ^Jml + {pp b + p^-^/2. The high energy data 
from UAl and CDF [35| are also shown in Fig. |2] and agree very 
well with the energy dependence of the color evaporation model, as 
obtained from eq. (§). Both sets of PDFs give compatible predictions 
at the LHC energy. The predicted cross section is up to a factor of 20 
larger at 5.5 TeV than that given by the empirical parameterization, 
da/dy\ y=0 oc exp(-14.7M T /v^) [g|, 0, labeled CR in Fig. @. 

The second step is to extract the T, T' and T" components separately 
using the T'/T and Y"/T ratios quoted earlier, and the known branch- 
ing ratios B^ i = S(T, — > where Tj represents the individual T 
states. If we write darjdy\ y=0 = fr t da / 'dy\ y= o for the Tj cross sections, 
then from eq. (0) 

f r B T + f r ,B r , + f r ,,B r » = B = 0.0016 . (10) 

Using fy/fx = 0.53 and fr"/fr = 0.17, as quoted earlier, and B^i 
from the Particle Data Book |4]§, we find 



fr = 0.044 , f r , = 0.023 , / T « = 0.0074 . (11) 



4 At high xf, additional production from other mechanisms such as intrinsic heavy 
quarks may be important. 
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10 L 



10" 




MR$ D- 
GRV HO 



10 J 



10' 



10^ 



10' 



Vs (GeV) 



Figure 2: We show Bda/dy\ y= o for x+T'+T" in pp collisions, as indicated 
in eq. (|), for the MRS D-' § (solid) and GRV HO @ (dashed) parton 
densities. Also shown is the empirical Craigie parameterization jl^, 42] 
(dotted). The data are taken from |53|, 34, 36, 38, 39, 41]. The vertical 



solid line indicates the maximum LHC Pb+Pb energy. This figure was first 
presented in Ref. p7y. 



• Thirdly, we must separate direct from indirect production. The mea- 
sured Tj production cross sections, or equivalently the fa., are only 
effective values which reflect both direct production and chain decays 
of higher mass states. Isolation of the direct production cross sec- 
tions for each Tj requires detection of the radiated photons associated 
with chain decays. Data which allows this separation is not currently 
available. As discussed later, the extent to which the required photon 
detection will be possible at the LHC is also uncertain. We shall see 
that this is a very important issue. For now, we make the assumptions: 

fxbi(iP) = /t' ; f Xbi (2P) = h" > ( 12 ) 

for each i, where i = 0, 1, 2 labels the Xb states in PDG notation ||44| ], 
and the d superscript indicates the / for direct production. With this 
and the summed branching ratios [Q 

J2 £?0fe(LP)^r 7 )~0.63, 

i=0,l,2 
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£ P( Xbi (2P)^T 7 )~0.16, 

i=0,l,2 

X; B(x w (2P)-T' 7 )~0.42 : 



i=0,l,2 



along with 5(T' -> TX) ~ 0.27 and T" decays to T', T and Xw(2P), 
we can compute the f d 's for direct production of all the Tj. One finds 
that about 0.014 of the total fy = 0.044 is due to xm(1P) decays. This 
would be similar to what is found in the analogous ip, \ c situation. 
Similarly, ~ 0.001 of fy and ~ 0.003 of fy> would be due to x fci (2P) 
decays. Also important are T' — > TX decays; fy [eq. (|TT])] implies that 
an additional 0.006 of fy would be indirect. The contributions from 
chain decays of the T" are small. Also, there are no contributions to 
the T" rate coming from higher states that are known to be significant. 
Altogether, we would have: 



f4 ~ 0.023 , ~ 0.020 , f4„ ~ 0.0074 
f x rf M(1P) ~ 0.023, f Xu{2P) 



/i(iP)~ 0.023, /' ~ 0.0074, (13) 



where i = 0,1,2 labels the different x^(lP, 2P) states. Note that only 
about half of fy is due to direct T production. The results in eq. fljjD 
imply that about 14% of a goes into bottonium state production. 

In Table f|, we give the corresponding normalized direct production cross 
sections in pp collisions, f d a, for each state with f d from eq. ( |T3"D and a 
computed using the MRS D— ' parton densities at 5.5 TeV. The cross sections 
given in Table |] are integrated over all rapidity. At this energy, the results 
with the GRV HO distributions are about 20% smaller. 

We would like to comment on the differences between the color evapora- 
tion model used here and other models of quarkonium production. The color 
singlet model f45f , in which bound states with the right quantum numbers are 



directly produced, has been shown to be insufficient to reproduce the Teva- 
tron data on J/ip, ip', and T production [|9|, [36|| . The low-pr Y production 
data from the Tevatron has recently been explained using a model of color 
octet production involving corrections to the bound state wavefunctions at 
higher order in the relative velocity of the QQ where a s ~ v for the bound 
state f4"6| . We would like to point out that the p^-dependence of quarkonium 



production is quite similar to the pr-dependence of unbound QQ pairs ]26). 



This same dependence also appears in the color octet model since the octet 
mechanism makes other 3 Si production channels available. Thus the color 
octet model is quite similar in spirit to the color evaporation model consid- 
ered here. However, at larger p?, p? > my, fragmentation processes such 
as g — ► T and b — > T j|7j should eventually dominate the production pro- 



cesses considered in the color octet model. It has been recently shown that 
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T 




y// 


Xb(lP) 


Xb(2P) 


f d a (nb) 


426 


233 


137 


426 


137 


iV d (PbPb) 


2.02 x 10 6 


1.75 x 10 6 


6.47 x 10 5 


2.02 x 10 6 


6.47 x 10 5 


iV^(PbPb) 


4.99 x 10 4 


2.30 x 10 4 


1.17 x 10 4 







Table 4: The normalized cross sections, f d a, for directly produced bottonium 
states in pp collisions at 5.5 TeV, using the f d 's of eq. (|T3f) and the prediction 
of eq. @ that a = 17.35 fib ^Z7§, where the MRS D— ' parton densities are 
employed with m& = 4.75 GeV and ft = rn T ^. Also given is the number N d 
of each type of bottonium state produced directly in Pb+Pb collisions at the 
LHC, calculated using eq. ([14]) with <r NN = 60 mb, T PbPb (0) = 32 mb -1 and 
L^ Pb = 2.59 nb -1 . For T, T', T" we give the corresponding number of fi + fT 



pairs from decay of directly produced states. 



the Tevatron charmonium and bottonium data is also in good agreement 
with the color evaporation model at NLO, as described in P5[ . In the color 
evaporation picture, gg scattering followed by the splitting g — > bb incorpo- 
rated at NLO is similar to models of g — *■ T fragmentation K7||. By including 
this splitting, the color evaporation model provides a good description of the 
quarkonium px distributions. 

We can now employ the results for f d given in Table ^ to predict rates 
for direct production of bottonium states in Pb+Pb collisions at the LHC. 
For central collisions, the expected rates are given by 



N d = a NN T PhPh (0)f d dL\ 



hit 



(14) 



where 0"jvjvTpbPb(O) is the number of central Pb+Pb collisions^ and f d a, with 
f d taken from eq. fljjD , is the cross section for direct production of a given 
bottonium state in pp collisions. In one month of running the integrated 
luminosity for lead beams is expected to be i^t Pb = 2.59/nb. Typical rates 
are given in Table they are on the order of 10 6 for T and T'. Approximately 
10-15% of the cross section is within \rj\ < 1. The number of muon pairs from 
the T, T' and T" decays, found by multiplying the total number of T, directly 
produced, N d , by the appropriate branching ratio, is also given in Table |4j. 
These rates suggest that production and suppression of these states should 
be measurable by CMS in the very clean fi + fi~ final state decay mode. The 
rapidity distributions are rather broad and nearly constant out to y ~ 4 with 
the MRS D— ' parton distributions. The GRV HO results exhibit a somewhat 
narrower rapidity plateau p7 |. 



5 We use <jnn ~ 60 mb, assuming that the high energy rise in a P p will also hold for 
crj\rjV) and IpbPb(O) — 32/mb. 
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The rates of Table ^ are those obtained by simply multiplying the pp cross 
section by the expected number of central Pb+Pb collisions. This procedure 
does not include any nuclear effects (e.g. shadowing or absorption) on T 
production. However, significant nuclear effects have been observed in fixed- 
target interactions. The E772 collaboration has measured a less than linear 
A dependence for T production at y^s = 38.8 GeV, aJ A = crJ p A a where 
a = 0.96 P5| . If no nuclear effects were present, the production cross section 
would grow as A 1 . Note that a = 0.96 implies a 20% reduction in the 
pPb production cross section relative to the A 1 prediction. The two possible 
sources of this reduction in pA collisions are shadowing and absorption. 

• If the reduction is due to shadowing, the observed greater suppression 
(smaller a) of J/ip compared to T implies that the shadowing of the 
gluon PDFs is greater at smaller x values. This would then imply that, 
for the T (or any other given onium resonance), shadowing would be 
greater at higher energies due to the lower x region probed. In heavy- 
ion collisions at LHC energies the shadowing mechanism could reduce 
crJ bPb by 25-40% relative to the result obtained assuming proportion- 
ality to A 2 ||49|| . The larger shadowing effect expected for Pb+Pb col- 



lisions compared to pPh collisions comes from the convolution of two 
nuclei with modified parton densities. To the extent that shadowing 
in central collisions is the same as for minimum bias collisions (which 
include peripheral as well as central collisions), the same 25-40% reduc- 
tion should be applied to the numbers of Table |J However, nucleons 
near the nuclear surface might experience a smaller shadowing effect 
than those in the center of the nucleus. 

On the other hand, the deviation of the A dependence from unity in 
fixed-target interactions has also been attributed to absorption by nu- 
cleons; interactions with comoving secondaries are possible although 
rare in pA collisions. In this approach, a is related to the the ef- 
fective onium absorption cross section in nuclear matter by 1 — a ~ 
9c a bs/167rrQ, assuming a step function nuclear shape. The apparent 
T absorption cross section is 40-50% less than the ip absorption cross 
section, suggesting that the T suffers fewer final-state interactions than 
the J/ij) (possibly because of its smaller size). The onium- nucleon ab- 
sorption cross section has been calculated on the basis of the operator- 
product expansion [j50[ ^TJ and has been shown to quickly attain its 



asymptotic value [52]. Absorption by nucleons has been shown to be 



less effective in peripheral relative to central nuclear collisions, see e.g. 



31| 20]. Therefore nucleon interactions with the J/ip may be insuf- 
ficient to account for J/ip measurements in heavy- ion interactions at 
current energies |[R [L6|. Secondary particle production increases in 
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central relative to peripheral collisions, increasing the influence of co- 
movers in central collisions as well as in nuclear collisions compared 
to pA interactions . The onium-comover cross section can be esti- 
mated from the onium-nucleon absorption cross section by the additive 
quark model assuming the comovers are mesons p3| . 



A mixture of shadowing and absorption is most likely. If the value of a 
depends on the Upsilon state, absorption is important: cr a b s could increase 
with the size of the Upsilon state ||54j| . If shadowing were dominant, a would 



remain essentially unchanged since the x values probed by the different Up- 
silon states are very nearly equal. In fact, the E772 collaboration observed 
no difference in the J/if> and if)' A dependence or in the T and T' A depen- 
dence within their uncertainties |48| . Additionally the if}'/ if) ratio appears to 



be independent of A |^7| but the data is such that a size-dependence of the 
absorption cross section cannot be ruled out. At the E772 energy, yfs = 38.8 
GeV, the target momentum fractions probed are x < 0.08 for the charmo- 
nium states, just in the shadowing region, while for the T states, x < 0.25, in 
the EMC or antishadowing regions. At the lower energies of other J/ if> mea- 
surements, shadowing should be less effective than absorption ||55|| . The mass 
scale in the PDF is also important if the shadowing is Q 2 dependent ||56|| . In 



Q 2 -dependent shadowing models, the gluon distribution evolves fastest with 
Q 2 ; the valence and sea quark evolution is rather weak for 2 < Q 2 < 10 
GeV 2 |49], [56]]. We note that nuclear deep-inelastic scattering data do not 



exhibit strong Q 2 dependence for Q 2 < 100 GeV 2 |57[| , in agreement with 
these models, but it is more difficult to probe the Q 2 dependence of gluon 
shadowing. Precision pA data is needed to sort out these effects, particu- 
larly at higher than current energies since at RHIC and LHC we will be in 
a new regime where x is small and Q 2 is relatively large. Measurements of 
dilepton production will prove helpful, as we discuss later, since the dilepton 
continuum could be influenced by shadowing but not by absorption. 

The pr-dependence of the suppression from shadowing/ absorption is es- 
pecially crucial for out later considerations. Since pA studies are not possible 
at the LHC, the exact rapidity and p^-dependence of shadowing/absorption 
suppression at the low x values probed at 5.5 TeV will be imprecisely known. 
Thus, it is fortunate that the pr-dependence of such suppression should be 
much weaker than that from QGP suppression (as we shall discuss below). 
As noted above, if shadowing is dominant, then suppression from this purely 
nuclear source should also be essentially independent of the Upsilon state. 
As our ensuing discussion will show, this would be a very important advan- 
tage. In contrast, if absorption is dominant the magnitude of the suppression 
could depend on the Upsilon state. Finally, we note that any broadening of 
the pr-dependence due to initial-state interactions in the nuclear target JSBfl 
may be expected to be negligible. Thus, we believe that any alteration of the 
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PT-spectra for the various members of the T family due to QGP effects will 
stand out clearly. Even more crucially, in the absence of QGP suppression the 
ratios of rates for different T family members should exhibit px dependence 
that is precisely that predicted without including shadowing or absorption.^ 

The pt- dependence for such ratios is then a direct probe of the QGP physics. 



4 Sensitivity of quarkonium suppression to 
the QGP model 

As we have just shown, the T production rate is high enough at the 
LHC for any QGP effects on the T family to be measurable. Our discussion 
in Sec. 2 suggested that at the high temperatures predicted for the minijet 
plasma, T production may be suppressed at the LHC. We would like to use 
this suppression to determine the initial conditions of the QGP, including 
T , t , and the temperature dependence of the QGP screening mass p(T). 
A simple measurement of the overall suppression of the production rates 
for the Upsilon family members will not be adequate. We must consider 
additional observables that could be affected by the presence of the QGP. In 
this section, we will show that the p^-dependence of the suppression provides 
a remarkable amount of information. 

If the QGP were of infinite temporal and spatial extent, resonance pro- 
duction would be suppressed as long as T > T D . But, there would be no 
significant px- dependence. In fact, however, both the lifetime and the spa- 
tial extent of the QGP are finite. Since a color singlet QQ pair takes some 
time to form a bound state, the competition between the formation time, 
the time needed for the QQ to separate to the bound state radiusf], and the 
characteristics of the plasma all affect the pr-dependence of the T states. 
We will now discuss the sensitivity of the onium suppression patterns as a 
function of pt to these ingredients. 

The finite lifetime of the plasma sets an upper limit on the pt at which 
the onium states can be suppressed. In the plasma rest frame, the bb forms a 



bound state in the time t F = ^r F where 7 = J 1 + (px/M) 2 (M is the onium 
mass) and the formation times, 7>, for the bottonium and charmonium states 
were given in Table ^| When the plasma has cooled below Tp, the production 
of the bound state is no longer suppressed. The time at which this occurs can 
be determined from the Bjorken model [59[] which assumes that the evolution 



of the system is described by ideal hydrodynamics, cooling by an isentropic 



6 If absorption is significant and <7 a b s depends upon the Upsilon state, then the ratio of 
PT-integrated cross sections would be affected but not the px variation. 

7 If the QQ is produced as a color octet, it is less likely to produce quarkonium than 
heavy- flavored hadrons. 
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expansion so that 



(15) 



Since the entropy density, s, is proportional to T 3 , the time at which the 
temperature drops below Tp is 



As long as io/tp > 1, quarkonium formation will be suppressed. Thus the 
maximum px at which suppression occurs is 



The values for to and pr.max for all the onium states are given in Tables [5] 
and ^| for both forms of /i(T), eqs. (|5|) and (^), and three different sets of 
initial conditions for LHC and RHIC energies. 

It is easy to check whether or not the bottonium states are suppressed 
given the formation time and T from the minijet initial conditions. The 
SU(iV) parameterization of eq. (||), with T c = 260 MeV, results in t D /r F > 
1 for all the members of the Upsilon family. Thus, all the states will be 
suppressed up to rather large values of px, suggesting that any suppression 
should be easily observable within the pt acceptance of CMS. However, To 
is reduced to 900 MeV when quark degrees of freedom are added as in the 3- 
Flavor scenario of eq. (||); the result is that tx>/ Tp < 1 for all the bottonium 
states and no suppression occurs. In other words, if the 3-Flavor quark- 
gluon plasma is finite and quickly equilibrated, none of the T states will be 
suppressed, and the p^-dependence should be unchanged from pp production 
modulo any nuclear effects on T production. 

In contrast, the charmonium states are suppressed in both the SU(iV) and 
the 3-Flavor /i(T) scenarios. The SU(iV) case gives values of px,m&x similar 
to or smaller than those of the bottonium family because of the smaller 
charmonium masses. This should be contrasted with the 3-Flavor QGP for 
which the J / ip would be suppressed over a very narrow range of pt relative 
to the Xc and ip' since the value of Tp for the J/ ip is nearly twice that of the 
other charmonium states (see Table [3] for To). If the px-dependence of the 
nuclear effects on the charmonium states is weak, then the ratio of higher 
charmonium states to the J/if> as a function of px would still be useful. 

We discuss the other initial conditions shown in Tables [| and [] and their 
consequences later. 

In a finite system, such as those produced in heavy-ion collisions, the 
suppression also depends on the size of the system. We present a schematic 




(16) 




(17) 
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LHC 




SU(JV) only, eq. (|5) 


gluons and quarks, eq. (6|) 


minijet initial conditions 




T = 1.14 GeV, t = 0.1 fm 


T = 900 MeV, t = 0.1 fm 




to (fm) 


PT,max (GeV) 


to (fm) 


PT.max (GeV) 


T 


2.49 


29.5 


0.105 







8.474 


43.5 


1.67 





Xb 


8.474 


30.7 


1.46 





V> 


8.474 


29.0 


1.54 


4.34 


ip' 


8.474 


20.6 


15.7 


38.6 


Xc 


8.474 


14.4 


17.5 


30.5 


Uncertainty relation initial conditions 


C = 1 




T = 600 MeV, to = 0.11 fm 


T = 785 MeV, t = 0.084 fm 


C = 3 




T = 346 MeV, t = 0.57 fm 


T = 453 MeV, t = 0.436 fm 




to (fm) 


PT,max (GeV) 


to (fm) 


PT,max (GeV) 


T 


0.39 





0.03 







1.3 





0.455 





Xb 


1.3 





0.82 





1> 


1.3 


3.27 


0.42 







1.3 





4.28 


9.9 


Xc 


1.3 





4.78 


7.6 


parton gas initial conditions 




T = 820 MeV, t = 0.5 fm 


T = 820 MeV, t = 0.5 fm 




to (fm) 


J"T,max (GeV) 


to (fm) 


PT.max (GeV) 


T 


4.6 


56.53 


0.19 







15.69 


81.98 


3.09 


12.83 




15.69 


58.9 


5.53 


18.58 




15.69 


54.0 


2.85 


9.34 




15.69 


38.5 


29.1 


71.6 


Xc 


15.69 


27.2 


32.5 


56.76 



Table 5: LHC values of to, and pT.max for the two choices, eq. (|5[) and eq. 
of fi{T), and for three sets of initial conditions. 
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RHIC 




SU(JV) only, eq. (|5) 


gluons and quarks, eq. (|6) 


minijet initial conditions 




T = 445 MeV, t = 0.1 fm 


T = 360 MeV, to = 0.1 fm 




*D (fm) 


£"r,ma X (GeV) 


t D (fm) 


£"T,max (GeV) 


T 


0.15 





0.003 







0.5 





0.05 





Xb 


0.5 





0.094 







0.5 





0.048 







0.5 





0.493 





Xc 


0.5 





0.55 





Uncertainty relation initial conditions 


C = 1 




T = 400 MeV, t = 0.164 fm 


T = 524 MeV, t = 0.126 fm 


C = 3 




T = 231 MeV, t = 0.85 fm 


T = 303 MeV, t = 0.65 fm 




t D (fm) 


PT,max (GeV) 


t D (fm) 


PT,max (GeV) 


T 


0.17 





0.013 







0.59 





0.2 





Xb 


0.59 





0.364 





V> 


0.59 





0.188 







0.59 





1.92 


2.94 


Xc 


0.59 





2.14 


1.31 


parton gas initial conditions 




T = 550 MeV, t = 0.7 fm 


T = 550 MeV, t = 0.7 fm 




t D (fm) 


J"T,ma x (GeV) 


t D (fm) 


£"T,max (GeV) 


T 


1.95 


22.3 


0.082 







6.63 


33.4 


1.3 





Xb 


6.63 


23.2 


2.33 







6.63 


22.6 


1.2 


2.8 


$ 


6.63 


15.9 


12.3 


30.12 


Xc 


6.63 


11.1 


13.7 


23.75 



Table 6: RHIC values of to, and pr,max for the two choices, eq. (0) and eq 
of fi{T), and for three sets of initial conditions. 
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calculation of the p^-dependence of the suppression following Ref. 50 . Sim- 



ilar analyses of QGP production have also been done E2, 61, 62fl. We can 



assume that the initial entropy profile is dependent on the radius of the QGP, 





so that 

tD(r)=t D (0)(l-\j i ) j , (19) 

where £d(0) is the value of tp> calculated in eq. ([16]) for resonances produced 
in the center of the system and R is the transverse radius of the plasma, 
assumed to be no larger than the radius of the nucleus, R ~ 1.2 A 1 / 3 . At any 
time t during the evolution of the QGP, the spatial boundary of the screening 
region is located at the transverse radius rs given by tp,(rs) = t, 

^^-fcwrf 2 ■ (2o> 

In general, if the QQ pair is created with transverse position x M = (0, r, 0) 



and momentum p^ = (J M 2 + py,p T ,0), at the formation time Tp the pair 



will form a bound state at x** = (tf\J 1 + (pt/M) 2 , r + T F p T /M, 0). If |r + 
t fPt/M\ > r s , the QQ pair escapes the QGP before t f and will successfully 
form bound states. This inequality will be satisfied (and the pairs will escape) 
for a range of angles between r and px — < 6 < 9 max (r,pT) — provided pr 
is such that (r + r F p T /M) > r s . Here, # max also depends on Tp, M, and r s as 
detailed below. The p^-dependent probability that the QQ pair survives and 
is able to form an onium state, S(pr), is the ratio of the number of bound 
states produced by those pairs that escape the plasma to the maximum 
possible number of onium bound states that would be formed at the given 
Pt in the absence of the QGP: 

s(pT) = /.*W<r,p r) 
7r J drrpyr) 

where we parameterize p(r) = (1 — (r/R) 2 ) a and 

f 7T Z <-l 

6(r,p T ) = I cos" 1 z \z\<l , (22) 
z > 1 
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where 




2 _ I T F p T 



(23) 



z = 



2r F r% 



When z > 1 the pair cannot escape, whereas when z < — 1 the pair always 
escapes. 

The results for the survival probability S(pt), with a = 1/2 and /3 = 1/4, 
are shown in Figs. 0(a), 0(b), and [5](c) for the T D values implied by eq. ([5]) (see 
Table ||) and for a decreasing sequence of radii: (a) R = i?pb = 7.1 fm, (b) 
R = R$ = 3.8 fm and (c) R — 1 fm, respectively. In all three cases, we observe 
that S(pt) ~ 0, corresponding to essentially total suppression, at small pr 
with our assumption of an equilibrated system. For R = i?pb the system 
is large enough for the finite lifetime to be more important than its finite 
size. Therefore S(px) = 1 for p T > pr,max while states with p T < pr,max are 
suppressed. If the QGP is of smaller spatial extent, the survival probability 
can become unity for px < PT,max since a smaller pt is sufficient for the pair 
to escape the system. This is illustrated by the results for R = R$ and R = 1 
fm displayed in (b) and (c) of Fig. |3|. As R decreases, the largest px value 
for which suppression occurs decreases more rapidly for the T' and Xb states 
than for the T because of their larger breakup time, tp] see Table [|. In fact, 
when R = 1 fm, S(pt) = 1 for the T' at pt ~ 10 GeV while the T remains 
suppressed up to pr ~ 20 GeV, as shown in Fig. 0(c). 

We should compare these results to other possible QGP scenarios. 

• First, there is the possibility that we have overestimated the amount 
of screening that occurs for a given level of minijet production. As we 
have noted, see Table |5], the screening times, to, are generally shorter 
for the perturbative estimate of eq. @ [even if quarks are not included 
by setting rif = 0] than for the SU(iV) gauge parameterization given in 
eq. (U)0. The result is that there is no suppression of the T states for 
the 3-Flavor QGP in equilibrium if minijet initial conditions apply. 

• It is also possible that the minijet density is different than we assumed 
in our calculations. For example, if the appropriate minijet momentum 
scale, po, were larger (smaller) than the value of po = 2 GeV assumed 
above, t and, consequently, the minijet density, T and, from eq. (|16|). 
t D would all be smaller (larger). A lower minijet density could produce 
a different suppression pattern for the T family; e.g. a choice of po 
is possible such that T production is not suppressed, while T' and Xb 
production are. 

8 The tp' and \c are an exception because the lower value of T c found when quarks are 
included reduced the breakup temperature, resulting in a longer screening time. 
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Figure 3: For the SU(iV) /x(T) of eq. (||) and the associated minijet initial 
conditions of To = 1.14 GeV, to = 0.1 fm, the suppression factor S(pt), 
calculated according to eq. (f2l|), is shown for T (solid), T' (dashed), and 
Xb (dot-dashed) production in three cases: (a) R = Rpb', (b) R = Rs; and 
(c) R = 1 fm. 



We note that typical early estimates of the initial conditions without 
minijet production, to = 1 fm and T = 250 — 300 MeV at the LHC (see 
e.g. fl2|l), produce no bottonium suppression regardless of the form of 



MT). 

An alternative set of initial conditions is obtained using the uncertainty 



principle [|63| , 64 1. In this case, to(E) ~ C%c, with (E) ~ 3T for a 
thermal average. In previous calculations, C = 1 and 3 were chosen. 
Those studies, using T c = 200 MeV and assuming a 3-Flavor QGP, 
obtained T = 0.46 GeV and t = 0.42 fm for C = 3 and T = 0.81 
GeV and t = 0.08 fm for C = 1 @]. In Tables || and || we update 
these results for the pure gluon SU(iV) and 3-Flavor QGP cases, using 
the appropriate T c for each. As can be seen from the formulae of 
Ref. H64R , the 3-Flavor scenario has a higher T than the SU(iV) case 
given our assumption that the final multiplicity is the same regardless 
of the initial conditions. We note that either choice of C produces 
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the same results because T oc yS/C (as shown in f6H| ) and t T oc 
C/3, implying that the C dependence cancels in the expression for to, 
eq. ([16]). Neither scenario leads to any bottonium suppression, and 
suppression of the charmonium states is only possible for rather low Pt 
values. 

The final initial condition scenario that we discuss is based on a parton 
gas model as adopted for J/ip suppression at RHIC and LHC energies 
in Ref. [p5|| . In this model one finds that the kinetic equilibration 
time is longer than for any of the cases discussed so far, t ~ 0.5 — 
0.7 fm. This time is reached when the momentum distributions are 
locally isotropic due to elastic scatterings and the expansion of the 
system. Chemical equilibrium is generally not assumed but the system 
moves toward equilibrium as a function of time. In this picture, the 
cooling of the plasma is more rapid than the simple Bjorken scaling 



picture []59| we have adopted, producing incomplete suppression at low 
Pt- To make a schematic comparison with our assumption of complete 
equilibrium, we give in Tables || and |6] the values of t D and pr.max for the 
onium states computed using the same T and t as employed in their 
calculation, derived from the HIJING Monte Carlo code p6[ . (Note 



that we use the same initial values independent of the form assumed 
for /i(T).) Because the equilibration time of the parton gas is longer 
than that obtained from the minijet initial conditions, the time the 
system spends above the breakup temperature is also longer, leading 
to suppression for both the SU(iV) and the 3-Flavor forms of /x(T), 
even though To is lower. These results are illustrated in Fig. |] for the 
two /u(T) parameterizations. Figs. |4](a) and|](b) illustrate the survival 
probability S(px) for the SU(iV) plasma assuming R = R Ph and 1 fm, 
respectively. In the SU(iV) case, pr,max is sufficiently large that the lead 
nucleus is already small enough for the pt range of the suppression to 
be less than the calculated pr.m&x- The 3-Flavor QGP results are shown 
in Figs. ^|(c) and |](d) for R = Rp\> and 1 fm, respectively. Here the 
temperature is not high enough for T suppression but the T' and Xb 
are suppressed, albeit over a much smaller px range than in the SU(iV) 
case. 

A given prediction for the survival probability S(pt) can be used to cal- 
culate the pt distributions of the T states. Since S(pt) depends upon the T D 
and Tp of each Upsilon resonance, each state will have a unique pt depen- 
dence, despite the fact that the pt distributions determined from the color 
evaporation model are the same for directly produced T, T' and Xb before 
inclusion of QGP effects. As we have noted, the precise S(pt) factors for 
the different T resonances also depend on the choice of /x(T), to, T , and R. 
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Figure 4: For parton gas initial conditions of To = 820 MeV, to = 0.5 fm, 
the suppression factor S(pt), calculated according to eq. (|2l|) , is shown 
for T (solid), T' (dashed), and \b (dot-dashed) production in four cases: 
(a) R = R Ph and /x(T) of eq. (@); (b) R = 1 fm and /z(T) of eq. (|); (c) 
R = i?pb and n(T) of eq. (^); (d) R = 1 fm and n(T) of eq. @. 



The resulting spectra are shown in Fig. ^] in the case of S(pt) taken from 
Fig. |3|(a), assuming the SU(iV) form of /x(T), with R = i? Pb , t = 0.1 fm 
and To = 1.14 GeV. The upper curves are the pt distributions in pp colli- 
sions at y/s = 5.5 TeV while the lower ones are the pp curves multiplied by 
S(px) to indicate the pt distribution expected in the presence of the quark- 
gluon plasma. Note that plasma production decreases the slope of the pr 
distribution. 

The raw px distributions before including the survival probabilities were 
obtained using the color evaporation model predictions calculated with the 
next-to-leading order QQ code of Ref. J2BJ. It is well-known that at low 
Pt it is necessary to resum the next-to-leading ln(M 2 /p^) powers for each 
power of a s to get an exponential factor. This has been extensively stud- 
ied for W and Z production |)7|, [68|, ^9], [7C], |7T], 0. The procedure in- 



volves computing a form factor in impact parameter (b) space, and then 
Fourier transforming back to px space. For Aqcd *C 1/b < M, the form 
factor can be computed perturbatively; for b > 1/Aqcd non-perturbative 
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Figure 5: The pr distributions for (a) T, (b) T', and (c) Xb production 
are given for pp interactions (solid) and plasma production (dashed). The 
plasma results are the pp distributions multiplied by S(pr) as computed for 
R = Rpb, h = 0.1 fm and T = 1.14 GeV and assuming the SU(iV) form 
of n(T) as given in eq. (||); S(pr) for this case was plotted in Fig. |3](a). We 
have taken y/s = 5.5 TeV. The differential cross sections (before including 
S(px)) are normalized to the integrated values tabulated in Table |j. 



modeling is necessary. In the latter range, the conventional procedure is to 
introduce a cutoff in b and a non-perturbative factor in the fe-space form 
factor. It is found [72] that the non-perturbative component dominates the 
PT-dependence modifications for < 6 GeV, implying substantial uncer- 
tainties for ^-dependence in this domain. However, Ref. 



72 also demon- 



strates that the ratio [da(W) / dpT\/[d<j{Z) / dpx] is very insensitive to the 
non-perturbative uncertainties due to the fact that mw is not very different 
from rriz- Further, the ratio deviates only slightly from the perturbatively 
computed value. Applied to the closely-spaced resonances of the Upsilon 
family, the formalism of Ref. |72| implies that the color evaporation model 
predictions for the ratios of the production rates for different T states will also 
not be significantly modified by the non-perturbative physics. The ability to 
reliably compute a ratio such as T'/T as a function of pt down to relatively 
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low pt values is critical to the program described in the next section. 

5 A program for determining the nature of 
the QGP 

We have seen that the uncertainties in the initial conditions as reflected 
in to and To, the uncertainty in the QGP radius R, and the ambiguity in the 
temperature dependence of the screening mass, fi(T), all imply that there 
are many possible results for the relative suppression of the T states. The 
challenge will be to extract the correct plasma picture from the data. We 
describe below the inputs, procedures and detector features required in order 
that the production and suppression of the T states can be used as a tool to 
determine the true state of the plasma and the proper initial conditions. Data 
from the LHC and RHIC will prove complementary; it is likely that Upsilon 
family production rates must be measured at both machines in order to 
reach reliable conclusions. Further, measurements of the charmonium states 
at both LHC and RHIC will provide valuable information. We organize the 
following discussion into a series of (overlapping) tasks and requirements 
necessary to extract information about the QGP. 

Checking the color evaporation model and establishing a baseline. 

As a first step, the production rates and pr spectra for the T states must 
be measured in pp collisions at the standard LHC nucleon-nucleon collision 
energy of 5.5 TeV and, by lowering the energy of the LHC still further, at 1.8 
TeV.[] This will allow us to check the predictions^ of the color evaporation 
model. In particular, we must confirm that the the T'/T and T"/T ratios 
are independent of energy and pr and that the (common) pr-dependence of 
the individual cross sections agrees with the color evaporation prediction. 
These measurements will allow us to confirm our perturbative modeling and 
understanding of T production in the clean pp environment, thus establishing 
a baseline for comparison with results obtained for nuclear collisions. 

Determining the effects of shadowing and nucleon absorption using 
RHIC data. 

Next, we must determine if the T'/T ratio is independent of purely nu- 
clear matter effects (e.g. shadowing and absorption) by measuring the ratio 

9 The lower energy measurement should agree with that obtained at the Tevatron []36f 
to the extent that valence quark distribution function contributions do not contribute 
significantly in pp collisions, but a direct pp measurement will remove the necessity of 
making this assumption in order to carry out the following program. 

10 Assuming that gluon fusion is dominant at both LHC and Tevatron energies. 
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in pA collisions. Since these studies are not possible at the LHC, this must 
be established at RHIC. The heaviest nuclear beam at RHIC will be gold, 
but the shadowing should be quite similar in gold and lead. Since RHIC is 
a dedicated heavy-ion facility, the PHENIX experiment at RHIC is expected 
to obtain several thousand T events per year along with a similar number 



of T"s JF3]. This should be sufficient to determine: a) the px dependence 
of the T'/T ratio in pA interactions where nuclear effects are present but 
no QGP is expected; and b) the magnitude of the T'/T ratio. If the ratio 
is the same (smaller) than in pp collisions, then shadowing (absorption) is 
probably the dominant purely nuclear effect. The dependence of the onium 
ratios on shadowing/absorption may also be established by studying ip'/ijj in 
pA collisions: lower x partons are probed and the rates are higher. 

If shadowing is dominant over absorption, then measurements of the ac- 
tual magnitude of the gluon shadowing will be important if we are to use 
the absolute pr-dependence of the onium cross sections instead of their ra- 
tios to analyze the QGP. It has been suggested that pA measurements of 
the dilepton continuum (M > 1 GeV) at RHIC will provide a measure of 
gluon shadowing since dileptons from charm production and decay (mainly 
from gluon fusion) will be the largest contribution in the mass region below 
the J/ip [f74| . Such charm studies are particularly useful because they probe 
gluon shadowing at the low x values that will be important for Upsilon pro- 
duction at the LHC. Measurements at RHIC for dilepton masses near the T 
mass probe parton x values of x ~ 0.1 (at y ~ 0), i.e. substantially above 
those appropriate at the LHC and also in a range where little shadowing 
is expected In order to probe small x values for dilepton masses near 
Mr, it would be necessary to accumulate good statistics for the T's at high 
rapidity, which can be observed in the muon arms of the detector. An T with 
y ~ +2 or —2 will have £2 ~ 0.35, i.e. in the EMC region, or x 2 ~ 0.006, 
in the shadowing region, respectively. The latter would probe nuclear effects 
on T production for gluon x values that are dominant for T production at 
the LHC around y ~ 0. 

Perhaps the best way to determine the importance of nucleon absorption 
for the onium states is to measure the relative rapidities of the onium states 
and the initial nucleons since the baryon number in the central region of 
nuclear collisions at RHIC and LHC is expected to be small. If the nucleons 
and quarkonium states are far apart in phase space, the effects of nucleon ab- 
sorption should be negligible. A zero-degree calorimeter or a forward hadron 
spectrometer, such as that of the BRAHMS detector at RHIC [76 could 
clarify this matter. 

Determining the importance of comovers. 

The use of onium production data for determining the nature of the QGP 
will be difficult if there are non-QGP mechanisms (other than simple shadow- 
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ing) affecting onium production that are comparable to the predicted QGP 
effects. We have already discussed one such mechanism, namely nuclear ab- 
sorption. However, we concluded that it was relatively benign since it should 
have little pt dependence and will not be important if the projectile and tar- 
get fragmentation regions are well separated from the quarkonium states in 
the central region. A second mechanism of concern is the breakup/absorption 
of the onium states through interactions with comoving secondary hadrons. 
Indeed, the observed J ftp suppression [^, [L(| at yfs = 17GeV has been at- 



tributed to interactions with these comovers [20|. The comovers should also 
be present at higher energies. Such hadronic scattering can affect both the ip 
states and the T states. For example, the comovers can easily break up the 
J/ip if the combined comover-l- J/ip center-of-mass energy is above the DD 
threshold. Since the separation between the onium mass and the threshold 
is smaller for the higher onium states, the higher states will be more easily 
broken apart by comovers. Thus, comovers could cause ratios such as T'/T 
and ip' /ip to be smaller than expected. The T, being the most tightly bound 
of the onium states, should be least affected. Note that this is precisely the 
same as expected from nucleon absorption, which should be largest for states 
with large radii and small mass thresholds. As discussed earlier, if the co- 
movers are mesons, the onium-comover absorption cross section should be 
~ 2/3 the onium- nucleon absorption cross section. 

Although the interpretation of J/ip suppression at \/s = 17GeV as be- 
ing due to comover scattering implies high hadronic densities in the present 
experiments, a hadronic interpretation of quarkonium dissociation would be 
much less plausible at LHC energies, where, for all cases considered in Ta- 
ble [|, T is significantly above T c . The reasoning is as follows. For such high 
initial temperatures, it is almost certain that a QGP will form, delaying the 
creation of hadrons until the QQ pair or onium state has passed through 
much of the material. Thus, there would be little opportunity for interac- 
tions with comovers until the late stages of the collision. The To values for 
RHIC (Table |6|) are also generally higher than T c and significant breakup due 
to interactions with comovers again seems unlikely if a QGP is formed. 

While the above argument suggests that comover absorption will not sig- 
nificantly impact onium production, it will be important to try to determine 
the importance of comovers directly from experiment. We outline several 
possible approaches. 

• In pA collisions a QGP is not formed. Thus, as noted earlier, in pA 
collisions the dilepton spectrum in the M > 1 GeV region, being pri- 
marily due to cc production and decay, should provide a reliable probe 
of shadowing. Comovers will have only a small influence on open charm 
mesons because the scattering is elastic — a ttD interaction will change 
the momentum distribution of the D but will not take the charm out 
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of the system. Nucleoli-/} scattering is also expected to be elastic. On 
the other hand, a J / ^-comover interaction is likely to break up the J / ip 
into DD pairs. By comparing the 'shadowing' on and off the resonance, 
we can obtain a firmer handle on the magnitude of comover absorption 
of the J /ip. 

• In nucleus-nucleus collisions, the density of comovers which can absorb 
the onium state would be larger still, with the largest density in central 
collisions. Since these densities are related to the final multiplicity and 
thus to the transverse energy of the collision, perhaps one of the best 
ways to check the importance of comover breakup and absorption is to 
trigger on peripheral collisions (for which a QGP is not formed). One 
would compare predictions for the ip and T states based on shadowing 
in peripheral collisions with actual measurements to see if there is extra 
suppression due to comover breakup or absorption. Such collisions 
could provide very valuable guidance on whether comover absorption 
is important. In addition, one could compare (as above) the dilepton 
spectrum on and off the ip resonance. 

Comover absorption will also not occur if the onium px is large enough. 
Calculations of comover suppression as a function of px appear in Refs. W% 



|78| , |79| . The range of pt over which comover suppression occurs should not 
depend strongly on energy. One can make a crude estimate of the pt range 
of comover interactions as in eq. (17). In this case, 



PtVx = M v /(t co /r co )2 - 1 , (24) 

where t co depends on the size of the system, t co ~ R/v re \, (v Te \ is the relative 
velocity of the onium and the comovers, 0.6 for ?/>-comover interactions at 
the CERN SPS) and r co ~ 2 fm is the comover formation time in a nucleus- 
nucleus collision. Since we assume comovers are final-state hadrons, there is 
no To dependence and since the pt range depends on the comover formation 
time rather than the onium formation time, p™ max changes by 20% or less 
for each onium family. When R = Rpb, p™max = 5.8M. Note that if a QGP 
is formed, r co could be even longer since the comovers would not form until 
after the QGP hadronizes. At the very least, the above comparisons should 
allow us to determine the maximum px up to which comover absorption is 
important. 

The influence of comover absorption is clearly an important issue. In the 
following discussion, we assume that all nuclear and comovers effects on the 
onium production cross sections can be unfolded or proven to be small. 

Using RHIC data to establish the general nature of the QGP. 
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At RHIC ratios such as T'/T and tp'/ip will be measured in nucleus- 
nucleus collisions. For minijet initial conditions or initial conditions based on 
the uncertainty relations, the short initial time scales result in the prediction 
that there should be no significant suppression at RHIC for any of the onium 
states, regardless of the form of /x(T). However, in the parton gas model the 
long t implies that suppression of the T states is possible for the SU(iV) 
(but not the 3-Flavor) form of /i(T). Further, the ip states are predicted 
to be suppressed for both the SU(iV) and 3-Flavor forms of /x(T). Thus, 
the suppression pattern of the ip and T states will help distinguish the fJt(T) 
dependence if significant suppression is observed. Specifically: 

• If all the ratios are independent of p?, then either no plasma has been 
produced or the initial conditions of the plasma are unfavorable for 
suppression to occur. 

• If the ip'/ip ratio depends on px, but the T'/T ratio does not, the most 
consistent picture is that a QGP has formed with initial conditions like 
those of the parton gas model and /x(T) is of the 3-Flavor form given 
in eq. (|). 

• If both ratios are p^-dependent, then parton-gas-like initial conditions 
are required and the screening mass is more likely to be of the SU(iV) 
form as in eq. fl5|). 

Analyzing the QGP at the LHC. 

One advantage of the LHC over RHIC is that onium suppression can 
occur for a much broader range of possible initial conditions. In particular, 
at the LHC the minijet initial conditions produce an initial temperature T 
that is sufficiently large that the fast equilibration predicted (i.e. small to) 
does not prevent suppression. More generally, the higher T possible at the 
LHC will be crucial since fast equilibration is perhaps more plausible. 

The question we now address is how to best determine to and other fea- 
tures of the QGP. A high statistics study of the T'/T ratio as a function of 
p T may prove conclusive. If significant py-dependence of T'/T is found for 
Pt ft 10 GeV in Pb+Pb collisions at 5.5 TeV, then it will be virtually certain 
that a quark gluon plasma was formed. (As argued at the end of the previous 
section, any p^-dependence of this type of ratio arising from low-pr pertur- 
bative or non-perturbative resummation effects should be small, and effects 
from comovers/absorption/shadowing should be absent once pr is reasonably 
large.) The precise behavior of the T'/T ratio can then be used to strongly 
constrain the QGP model parameters. In particular, the ratio will be very 
different if only the T' is suppressed in comparison to the case where both 
states are suppressed. 
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Prompt Production Ratios 
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Figure 6: The direct or prompt T'/T ratio as a function of pt is shown 
for several choices of initial conditions at the LHC. In (a), we show results 
for the minijet model with the SU(iV) form of fJ>{T), eq. (0), (for which 
T = 1.14 GeV and to = 0.1 fm) taking R = R Ph (solid), R = R s (dashed), 
R = 1 fni (dot-dashed). In (b), we consider the parton gas model with 
T = 0.82 GeV and t = 0.5 fm, and show results for the SU(iV) fi(T) with 
R = Rpb (solid) and R = 1 fm and for the 3-Flavor fi(T) from eq. (g) with 
R = i?pb (dot-dashed) and R = 1 fm (dotted) . 



In Fig. H we show the T'/T ratio for several sets of initial conditions. 
The 'prompt' ratio displayed is that obtained if only directly produced T' 
and T states are included. In other words, the ratio is computed using /£, 
and fy from eq. ([13|) . We will return shortly to the issue of whether or not 
this prompt ratio can be experimentally measured. 

Fig. |6](a) displays results obtained when the minijet initial conditions are 
used with i — 0.1 fm and /x(T) is as given in eq. @. Both the T and 
T' states are suppressed and the solid, dashed and dot-dashed curves (for 
R = Rph, Rs and 1 fm, respectively) show that the pr-dependence of the 
ratio can distinguish different plasma radii. At low the T is actually 
more suppressed than the T', leading to a significantly larger ratio than that 
expected in low energy pp collisions (or found in current pp measurements). 
For a large R, the ratio eventually drops below the value expected in the 
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absence of QGP effects (i.e. for large enough p?). The point at which this 
shift occurs is at p T « 25 GeV if R = R Ph and at « 29 GeV for R = R s . 
The kink in these ratios at 30 GeV occurs at the point where the T is no 
longer suppressed. In a system with a smaller R = 1 fm radius, the T is 
always more suppressed than the T' so that the ratio returns to the pp value 
at p T « 20 GeV. 

Fig. |6|(b) shows the ratio computed using the parton gas initial conditions. 
When the SU(iV) fi(T) is assumed, the behavior is similar to the minijet case 
although when R = Rp^ the increase in T'/T persists to much larger px- 
When R — 1 fm, the result is almost indistinguishable from the minijet case 
with the same radius. The maximum of the ratio is ~ 10% larger at pt ~ 1—2 
GeV, but this may be difficult to distinguish. If the 3-Flavor [eq. @] scenario 
is assumed for the screening mass, then the T' is suppressed while the T is 
not and the pr-dependence of the ratio is quite different. The dot-dashed 
and dotted curves of Fig. ||(b) show that until the effects of the plasma 
disappear, the ratio is less than expected in the absence of QGP formation. 
However, it is more difficult to differentiate between a system the size of the 
lead nucleus or a system with R = 1 fm unless the measurement is made 
with high statistics. We have already noted that if significant suppression is 
observed at RHIC in the ip'/ip ratio as a function of pt and T suppression 
is observed at RHIC, the general behavior of fi(T) should be relatively clear 
by the time the LHC measurement is made. 

Finally, we recall that for the uncertainty relation initial conditions, both 
to and To are small even at the LHC and no suppression is expected for 
the Upsilon states. To learn about the existence and nature of a QGP via 
measurements of onium production will not be possible for such models. 

The importance of separating prompt from indirectly produced 
onium states. 

An important issue is whether or not the 'prompt' ratio discussed above 
can be extracted experimentally. We saw earlier that a significant fraction of 
T production arises from T' and Xb (Xb denotes the collective x^(lP, 2P), i = 
0, 1,2 states) production followed by T' — > TX and Xb lT, respectively. 
Further, some of the T"s come from xmC^P) - > jT'. We can only isolate 
the prompt ratio discussed above if we can veto such decay sources of the T 
and T'. In the case of the Xb decays, this will only be possible by detecting 
the associated photon and reconstructing the Xb- The ability of the LHC 
detectors to do this will be limited at low pr for two reasons: i) there is a 
large ir° background from which it is difficult to extract isolated 7's until 
the pt of the 7 becomes substantial (> 10 GeV) and ii) the photon energy 
resolution will be poor at low pt- In order to veto T' — ► TX decays we must 
be able to detect the X = 7r + 7r~, 7r°7r° decay products. This latter may be 
impossible; however, this indirect decay is only responsible for about 13% of 
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Production Ratios With Indirect Decays 
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Figure 7: The T'/T ratio including indirect decays as a function of px is 
shown for several choices of initial conditions at the LHC. In (a), we show 
results for the minijet model with the SU(iV) form of //(T), eq. (||), (for 
which T = 1.14 GeV and t = 0.1 fm) taking R = R Ph (solid), R = R s 
(dashed), R = 1 fm (dot-dashed). In (b), we consider the parton gas model 
with T = 0.82 GeV and t = 0.5 fm, and show results for the SU(iV) fi(T) 
with R = i?pb (solid) and R = 1 fm and for the 3-Flavor fx{T) from eq. (||) 
with R = i?pb (dot-dashed) and R = 1 fm (dotted). 



the T"s, so that if the Xb decays can be vetoed we will obtain a result much 
like the prompt ratio plotted in Fig. If the Xb decays cannot be vetoed, 
then we must consider all sources of the T' and T, each of which will be 
associated with a different suppression factor. Numerically, it is sufficient to 
consider: 

r + x »(2P)H t') + t"(-> v) 

T + Xb(l J P,2P)(^T) + T / (^T) + T"(^T)' 1 ' 

In computing this 'indirect' T'/T ratio we have assumed that the survival 
probability, S(pr), of the x&(2-P) states is the same as that for the Xb(lP) 
states and that S(pt) for the T" is the same as for the T'. 

The indirect ratio, denned in eq. ( p5[) is plotted in Fig. [F[ Clearly, the 
indirect ratio is somewhat less sensitive to changes in R and n(T) than is the 
direct T'/T ratio. If the SU(iV) form of /x(T) is appropriate, the figure indi- 
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cates that a measurement at the 20% level is needed to distinguish between 
the pp value of the ratio and the QGP prediction including indirect decays 
in both the minijet and parton gas initial condition scenarios. (If there are 
substantial systematic errors in the ratio, the detection of a deviation might 
be most difficult for the parton gas SU(iV) case when R = Rp^, since in this 
case the ratio varies slowly with p^-) If the slowly-growing 3-Flavor form of 
fJ>(T) is correct, the dot-dash and dotted curves of Fig. 0(b) show that the 
parton gas model predicts significant and strongly p r -dependent departures 
from the T'/T ratio expected in the absence of a QGP. It would even be 
possible to determine the size of the system using the indirect ratio — as 
shown in Fig. 0(b) the Xb contribution to the T causes the T'/T ratio to 
increase above the pp level at px ~ 12 GeV for R = Rp^ but not for R = 1 
fm. 

Thus, all is not lost if only the indirect ratio can be measured, but sensi- 
tivity is reduced and the accuracy with which the T'/T ratio can be measured 
becomes crucial. Also, our ability to predict theoretically or measure in other 
contexts with some reliability the relative values of the f d, s becomes quite 
important. 

Estimates of the statistics required and the importance of vetoing 
indirect decays. 

The ability to employ the T'/T ratio to analyze the nature of the QGP 
will depend crucially on the statistics available. For a NN Tp h p h (0)Lf^ Ph = 
4.973 x 10 3 nb _1 per month, as assumed in Table f|, and the da(pp) / dpx x 
S(pt) spectra of Fig. |[| which range very roughly from 1 nb/GeV to 0.1 
nb/GeV over the pt range of interest, we find ~ 5 x 10 3 falling to ~ 5 x 10 2 
events per GeV for the individual T, T' and Xb resonances. For branching 
ratios of £?t — 0.025 and Bx> = 0.013 and an efficiency/acceptance factor 
of 33% (which, as noted earlier, applies to the CMS detector), this leaves 
us with ~ 41.3 (~ 21.5) falling to ~ 4.13 (~ 2.15) T (T') events per GeV 
per month of running. For a 2 GeV bin size, 12 months of running and 
negligible background (see below), we obtain T (T') statistical accuracies of 
~ 3% (~ 4.5%) at low p T rising to ~ 10% (~ 14%) at p T ~ 25 - 30 for the 
individual spectra, and of 5.4% rising to 17% for the T'/T ratio. Thus, if 
the prompt ratio can be extracted, somewhat more than a year of running 
at Lf^ Ph = 2-59 nb _1 per month would be sufficient to measure the ratios 
with the required accuracy over the full range of px values of interest. If the 
T'/T ratio can only be measured including indirect decays, greater accuracy 
is required at large px- Roughly 36 months of running at Lf^ Pb = 2.59 nb _1 
per month would be required to achieve < 10% accuracy for the ratio at 
Pt ~ 30 GeV. Since this amount of running in the heavy ion mode at the 
LHC is quite unlikely, we see how important the vetoing of indirect decays 
is if we are to use the T'/T ratio to analyze the QGP. 
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The importance of mass resolution. 

Also crucial to the program is the ability to resolve the T, T' and T" 
resonances from one another. The CMS detector provides excellent mass 
resolution in the fi + fi~ mode, of order 40 MeV, adequate for clean separation 
of the T, T' and T" signals (see Fig. 12.31b in Ref. [|21|| ). As also shown 
in plj, the excellent resolution implies a very modest level of background 



under the T peaks. Preliminary estimates of the mass resolution for the 



ATLAS detector are less encouraging, of order 200 — 300 MeV pD |, for which 
separation of the resonances would not be possible. 

Other onium ratios. 

In the above discussion, we have focused on the T'/T ratio. Other ratios 
would provide complementary information. However, the T"/T ratio is a 
factor of 4 smaller [see eq. fll3"|) 1 at the direct, i.e. prompt, level and statis- 
tics would be a problem. The direct /prompt x&(lP)/T ratio is similar to 
the T'/T ratio, the major difference being that the x&(lP)/T ratio never de- 
creases below the pp level at high p T in the SU(iV) case since p^max ~ Prmax 
(see Figs. |3] and f|). Further, if the x&(lP) radiative decays can be detected 
with reasonable efficiency, the statistical accuracy with which the Xf,(lP)/T 
ratio can be measured would be comparable to that for the T'/T ratio, thus 
making it a most valuable addition. This adds to the importance of being 
able to trigger efficiently on the decay photons so as to isolate these two 
direct ratios. 

Extracting individual survival probabilities and using Z production 
to control systematics. 

So far, we have considered only ratios of onium rates. However, it will be 
important to know that suppression of the general type expected is taking 
place in the individual spectra contributing to these ratios. It would be even 
more valuable if a quantitative determination of S(pt) for the different reso- 
nances were possible through a direct comparison of the differences between 
the pp spectra and the Pb+Pb spectra at 5.5 TeV, as illustrated in Fig. |^. 
In principle, the ratio of the properly normalized Pb+Pb spectrum to the pp 
spectrum for a given resonance simply gives S^Pt^Q provided there are no 
non-QGP effects that can also alter the Pb+Pb spectrum. However, it is very 
likely that there will be some nuclear effects from shadowing, comovers and 
absorption. Effects from comovers and absorption cannot be disentangled 
except as described earlier. However, there is some possibility of obtaining 
additional information regarding shadowing. 

The strength of shadowing depends upon the x values at which the parton 
distributions are probed. Thus, the RHIC pA measurements of onia produc- 

n As argued in the previous section, non-perturbative resummation effects upon the pr 
distributions will cancel out in such a ratio. 
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tion, while helpful, cannot directly determine the full effect of shadowing at 
5.5 TeV where we wish to determine the presence of the QGP and measure 
S(pt)- An independent test of the shadowing effect using measurements at 
y/s = 5.5 TeV would be of great value for checking our understanding of 
shadowing and assessing the accuracy with which it can be unfolded from 
the data so as to reveal the QGP effects directly. 

pp^AiV"X at VeT=5.5 TeV 
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Figure 8: In (a), we plot dajdpT for pp — > /i + /i~ X at = 5.5 TeV as a 
function of pt for Z production (defined by mz — 5 GeV < M^+^- < mz + 
5GeV). In (b), the same cross section is plotted for 15 GeV < M^^- < 
20 GeV. The separate contributions from qq and qg + qg collisions are 
indicated by dashed and dotted histograms, respectively. 



In the current experiments at the CERN SPS, J/ip production is com- 



pared to the dilepton continuum fll^, |i6| . The continuum is assumed to be 



produced via the Drell-Yan (7*, Z*-exchange) process and is, in fact, Drell- 
Yan-like in its behavior. At RHIC energies and higher, the continuum will 
be more difficult to understand because of the important contribution from 
semileptonic cc and bb decays. Not only are there uncertainties in the to- 
tal cc cross section, but the heavy quark decays are also subject to nuclear 
effects. The relatively small Drell-Yan contribution is also subject to shad- 
owing effects in the mass range between the J/tp and the T. Another choice 
is needed. One possibility is Z production. Because the Z is produced in 
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point-like fashion, the difference between the pr-dependence of the Z in pp 
collisions vs. Pb+Pb collisions will not be influenced by the quark-gluon 
plasma^. 

Figure § illustrates the cross section for Z production as a function of 
Pt at = 5.5 TeV assuming no shadowing. There, we plot da/dpr for 
pp — > ZX and the separate contributions from qq and qg + qg collisions. 
(The distributions include 1-loop corrections to the Z + jet cross sectionP], 
but do not include \ow~pt resummation effects.) To obtain the number of 
events per month per GeV coming from central collisions^, we multiply by 
the factor cr NN T PhPh (0)Lf^ Ph = 4.973 x 10 3 nb" 1 . At p T = 50 GeV, the 
cross section is of order 10~ 3 nb/GeV, implying about 2 events per GeV for 
40% acceptance and detection efficiency. Thus, for a 5 GeV bin we obtain 
about 10 events in this bin per month of running. After a year of running, 
this would yield a statistical accuracy of order 9%. At low p T , event rates 
are a factor of > 10 larger, yielding correspondingly greater accuracy. The 
predicted effects of the QGP, see for example Figs. ^ and [|, typically imply 
survival probabilities that differ by much larger percentages compared to 
unity. In any case, as estimated earlier, the errors in the measurements of 
the Upsilon resonance spectra will be larger. Thus, production rates in the 
Pt ^ 50 GeV domain are high enough that Z — > e + e~ , fi + fi~ can provide a 
standard of comparison. Of course, we will also wish to measure Z production 
in pp collisions at \fs = 5.5 TeV so as to determine if shadowing and other 
nuclear effects have impacted the Pb+Pb spectra. Good statistical accuracy 
for up to pt = 50 GeV requires integrated luminosity for pp collisions of order 
L = 0.01 fb _1 , which should be easily achieved, even at this reduced energy, 
by running for a few weeks.0 

The two difficulties with using Z production as a benchmark are: (i) mz 
is much larger than the masses of the Upsilon resonances and (ii) the qq and 
qg + qg parton luminosities are probed vs. the gg luminosities responsible 
for bb production. 

• The large size of mz compared to the Upsilon family masses reduces 
the value of Z production as a benchmark for two reasons. First, shad- 
owing and related nuclear effects may well be dependent upon Q 2 . For 

12 Once again, this ratio will be very insensitive to low-pr perturbative and non- 
perturbative resummation modeling. 

13 We thank U. Baur for providing a program against which to check our own calculations. 

14 We wish to restrict observations to central collisions in case shadowing and nuclear 
effects (that might possibly be present for the Z) for peripheral collisions are somewhat 
different. 

15 If we assume that low-luminosity running yields pp integrated luminosity of 10 fb _ per 
year at y/s — 14 TeV and that instantaneous luminosity falls as E 2 , then the corresponding 
integrated luminosity at \fs = 5.5 TeV would be ~ 1.5 fb _1 per year or ~ 0.03 fb _1 per 
week. 



37 



example, in the limit of large Q 2 , holding the parton x fixed, shadowing 
becomes stronger in some models (see, e.g. , Ref. ]8l]] whereas in others 



it is predicted to disappear altogether (see, e.g. , Ref. |p2[|). See also 
Ref. p3| . Thus, it is very possible that the shadowing at Q 2 = M\ 



will differ substantially from that at Q 2 = m 2 z . Second, the x values 
probed (x ~ mz/ 'y/s ~ 0.016 at y = 0) are much larger than in Upsilon 
production at the same energy. (Detailed measurements as a function 
of the Z rapidity would provide information about smaller x values, 
but require more luminosity for good statistics.) 

The separate contributions of the qq vs. qg + qg densities to the Z 
spectrum appear in Fig. || We observe that qq collisions are dominant 
forpr ^ 15 GeV; for higher pt values qg+qg collisions dominate. Thus, 
to probe nuclear effects on the g distribution at Q 2 = m|, we must have 
knowledge of their effects on the q and q distributions for x values in 
the vicinity of those being probed at this same Q 2 . At low to moderate 
Q 2 the nuclear effects on q and q PDFs are rather well measured and 
do not appear to be strongly Q 2 dependent [^7]. However, there is 



currently no direct measurement of q and q shadowing at small x with 
Q 2 values as high as m 2 z . If nuclear beams become available at HERA, 
then direct measurement of shadowing at appropriate x and Q 2 values 
would be possible. 

Even though Z production does not directly probe gluon shadowing at the 
x and Q 2 values of interest for Upsilon production, Z production would be 
of considerable value in developing a reasonable theoretical understanding or 
model of shadowing. If we measure the pt dependence of Z production and 
find that it is consistent with one or more theoretically reasonable models of 
shadowing, this will already be a very important benchmark in that we can be 
confident there is not some totally unexpected physics occurring in nucleon- 
nucleon collisions. Of course, it is vital that Z production be measured both 
in Pb+Pb and in pp collisions at y/s = 5.5 TeV in order for us to develop 
confidence that we understand the impact of shadowing and other purely 
nuclear effects. 

Given the above issues regarding Z production, it would be advantageous 
if we could measure lepton pair production at pair masses nearer to those of 
the Upsilon states so as to constrain shadowing and nuclear effects at parton x 
and Q 2 values closer to those of direct relevance. We have already noted that 
we expect a large background from cc and bb production processes for lepton 
pair masses lower than the Upsilon family masses. At these low masses, this 
background will be very difficult to veto by requiring that the leptons be 
isolated because of the high density of soft tracks in the Pb+Pb collision 
environment. In the mass region above about 15 GeV the dilepton rate from 
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cc (bb) pair production is predicted to be smaller than (comparable to) that 
from 7*, Z*-exchange [Q. Further, in this higher mass range, vetoing the bb 
component using isolation requirements on the leptons might prove feasible 
at a level adequate to extract the pure DY dilepton spectrum. In Fig. [|, we 
have plotted da/dpx for production of muon pairs with 15 GeV < M^+ M - < 
20 GeV coming from 7*, Z*-exchange. We note that in this case the qg + qg 
collision component is always dominant, as would be desirable for learning 
as much about gluon shadowing as possible. However, the cross section is 
nearly a factor of 100 below that for production at the Z resonance, implying 
that statistics would be a factor of 10 worse. Even a year of running will not 
provide enough Pb+Pb luminosity to yield measurements that are sufficiently 
accurate to constrain the shadowing and nuclear effects at the needed level 
(< 5 — 10%). Thus, the low rate and uncertainty regarding our ability to 
veto the bb background imply that we should not count on being able to 
use dilepton pairs below the Z mass region to improve our understanding of 
nuclear effects on the gluon PDFs. Nonetheless, the possibility of doing so 
should not be ignored and appropriate data, including event characteristics 
that might allow vetoing, should be collected. 

HERA measurements of the nuclear g, q, and q densities as a function 
of x and Q 2 would greatly assist in a reliable extraction of S(pr)- Measure- 
ments overlapping the regions of interest for T and Z production at 5.5 TeV 
would be possible for an appropriate design of the facility. Using the HERA 
measurements for x ~ 0.02 and Q 2 = m|, it would be possible to unfold the 
shadowing effects in Z production and verify that the underlying px spec- 
trum is that expected on the basis of the Drell-Yan mechanism. This would 
allow a high degree of confidence that we could use the HERA measurements 
at the x and Q 2 values appropriate to onium production to unfold the un- 
derlying pt spectra of the onia as influenced only by absorption, comovers, 
and the quark-gluon plasma, with the latter hopefully having a much larger 
impact. 

6 Discussion and conclusions 

We have found that expectations for T (and also %fj) resonance production 
in Pb+Pb collisions at the LHC are very dependent upon the nature and 
details of the quark-gluon plasma. We have demonstrated that, far from 
being an impediment to our extracting the QGP physics, this dependence 
may very well allow us to determine much about the fundamental nature of 
the QGP, including: the energy density, the initial temperature, the plasma 
radius, and the temperature-dependent screening mass. Although we have 
kept the discussion of quarkonium production and suppression as general as 
possible, we chose the T family to examine in detail. The T rate is high 
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enough for statistically significant measurements to be made, particularly at 
the LHC. The T is also less affected by interactions with hadronic matter, 
providing a less ambiguous signal than the charmonium states. Strikingly 
different expectations for the pr-dependence of the T'/T ratio, the example 
upon which we focused, are found depending upon whether or not the QGP 
is formed and on the QGP properties. While our calculations of the survival 
probability are rather schematic, e.g. assuming full equilibration, they reflect 
the correct general trend for the scenarios discussed. 

The much higher energy of the LHC could be crucial for QGP suppression 
of onium production to be possible. This would be the case, for example, in 
the minijet model of the initial conditions. The high minijet density at the 
LHC implies that the pt spectra of the T resonances in this model are likely 
to be highly sensitive to details of QGP formation. In contrast, at RHIC, the 
minijet density will be a factor of 10 to 30 smaller, and the maximum initial 
temperature predicted is T ~ 0.5 GeV. Thus no suppression is expected. 
Initial conditions such as those predicted in parton gas models, for which 
the QGP takes longer to equilibrate, predict still stronger suppression at 
the LHC and could also lead to some suppression being observed at RHIC, 
depending upon the temperature dependence of the screening mass in the 
bb potential. There are also models of the initial conditions that predict no 
suppression even at the LHC. We point out that even a very high T does not 
guarantee plasma suppression. Models assuming a slow equilibration result 
in greater suppression because of the longer time the system spends in the 
screening region with T > Tq. 

Even if suppression is not present at RHIC, it is only at RHIC that pA 
measurements will be possible. Despite the lower energy, these pA mea- 
surements will greatly aid in constraining and checking our understanding 
of nuclear effects, such as those associated with shadowing, comover interac- 
tions and nuclear absorption. Thus, RHIC measurements will, as a minimum, 
provide important benchmarks. 

We have demonstrated that statistics should be adequate to detect the 
differences in the T'/T ratio that would discriminate between different QGP 
models (assuming that the ratio can be measured at the prompt level). Al- 
though this (and other) ratios have the advantage that various systematic 
effects, as well as the purely nuclear effects associated with shadowing, will 
cancel, further vital information can be extracted if the pr spectra for the 
individual resonances can be measured and the effects of ordinary nuclear 
shadowing on these spectra unfolded. We have pointed to a comparison with 
Z production as a potentially useful benchmark in this unfolding process. 
The ideal situation arises if, in addition to LHC Z data, HERA nuclear data 
on shadowing is available for the x and Q 2 values appropriate for both onium 
and Z production at the LHC. The observed p T spectrum for the Z could be 
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compared to the prediction based on the Drell-Yan process using the exper- 
imentally measured nuclear PDFs from HERA. If there is agreement, then 
we would have a high level of confidence that we can unfold the shadowing 
effects on the onia px spectra, using the PDFs determined from HERA data 
at the appropriate x and Q 2 values, so as to expose the impact of the QGP 
(which hopefully will be dominant over absorption and comovers). Whether 
or not the HERA data is available, Z production will provide an important 
benchmark as to how much nuclear shadowing effects alone are likely to be 
impacting the measured pt spectra of the onia. 

Certain detector features will play a key role in carrying out the analysis 
as envisioned here. Of greatest importance is the mass resolution in the fi + fi~ 
channel. The resolution for the CMS detector, as quoted in their technical 
proposal, is such that the T, T' and T" resonances can be cleanly separated 
from one another with very modest background under the peaks. Less certain 
is the extent to which the Xb - ► 7T process can be separated from direct T 
production by detecting the associated 7 for p^s in the range of interest. 
The statistics necessary to discriminate between different QGP models is 
much greater if the direct 'prompt' T'/T ratio cannot be isolated. One also 
would lose the ability to observe the Xb/"f ratio, which exhibits sensitivity 
to the QGP model similar to T'/T, thus roughly doubling the significance 
with which different QGP models can be differentiated. It will be important 
for the LHC detector groups to give further attention to both experimental 
issues. 
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